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I. INTRODUCTION 

The analyses of y-scaling (see, e.g. M) and super¬ 
scaling (based on the ^//-scaling variable) (see, e.g. [10l — 
[2l| l phenomena in inclusive electron scattering on nuclei 
have induced studies of (anti)neutrino-nucleus scattering 
on the same basis. This allows one to explore fundamen¬ 
tal questions of neutrino reactions and neutrino oscilla¬ 
tions in relation to the hypothesis of nonzero neutrino 
masses [22|. The theoretical concept of superscaling (a 
very weak dependence of the reduced cross section on the 
momentum transfer q at excitation energies below the 
quasielastic (QE) peak for large enough q and no depen¬ 
dence on the mass number) has been introduced yJJ, [HI 
within the relativistic Fermi gas (RFG) model. It has 
been pointed out in [li|, however, that the actual dy¬ 
namical reason of superscaling is more complex than that 
provided by the RFG model. This imposes the necessity 
to consider superscaling in the framework of theoreti¬ 


cal methods that go beyond the RFG model. An ex¬ 
ample is ea?., the Coherent Density Fluctuation Model 
(CDFM) [23|, [24| used in fl6l - fl9l [25j within this context. 

In [26| the analyses of superscaling have been extended 
to include not only QE processes but also those in which 
A-excitation dominates. The QE and A-region sca ling 
functions / CJE (V’ / ) and f A (ip r ) have been deduced in |26j 
from phenomenological fits to the data for electron-nuclei 
scattering cross sections by dividing the latter by appro¬ 
priate elementary N —»■ N and N —> A functions, respec¬ 
tively. Therefore they include all the effects of the nuclear 
dynamics, in particular NN correlations and final-state 
interactions (FSI), which should be reproduced by reli¬ 
able nuclear models. For instance, in 27 1 ]28| a QE scal¬ 
ing function with asymmetric shape has been obtained 
in agreement with the experimental scaling function us¬ 
ing the relativistic mean field (RMF) model for the final 
states. 

In order to exploit superscaling for neutrino-nucleus 
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studies, in [2f| (see also H3, H3) the above proce¬ 
dure has been inverted: the scaling functions have 
been multiplied by the elementary charged-current (CC) 
(anti)neutrino cross section to obtain the correspond¬ 
ing CC (anti)neutrino cross sections on nuclei for in¬ 
termediate to high energies in the same region of exci¬ 
tation. The scaling and superscaling ideas have been 
carried a step further in [13 to include neutral current 
(NC) (anti)neutrino scattering cross sections from 12 C, 
namely the reactions 12 C(iz,pzz)X, 12 C {D,pD)X involving 
proton knockout and 12 C(iz, nz^X, 12 C(zq nzz)X involving 
neutron knockout in the QE regime. The CDFM scal¬ 
ing function was applied to analyses of neutral current 
(anti)neutrino scattering on 12 C with energies of 1 GeV 
(u-channel inclusive processes) in [25| . A number of other 
theoretical studies have been devoted to both neutral- 
current (e.g. PH) and charged-current ( e.g. P-[4l|) 
neutrino scattering on nuclei. 

Our interest in the present work concerns the anal¬ 
ysis of different experimental data recently obtained on 
neutrino-nucleus processes in several facilities. At around 
1 GeV, data are available from the MiniBooNE collabo¬ 
ration, both for CC j42| and NCjEp neutrino- 12 C pro¬ 
cesses, and previous experiment [44j |. We note also that 
recent data for antineutrino-nucleus scattering are re¬ 
ported in phi , ITot . As known, the analyses of nuclear ef¬ 
fects in neutrino scattering are generally regarded as one 
of the main sources of systematic uncertainties in oscilla¬ 
tion experiments, in particular, when it concerns the un¬ 
derstanding of charged-current quasielastic (CCQE) in¬ 
teraction with nucleons bound in the nucleus {e.g. [471 ] 1 
in the energy range of around 1 GeV. For energies of 
a few GeV, the A-resonance excitation becomes equally 
important [481 ] (but keeping in mind that here the nuclear 
uncertainties are even larger). 

Though the main subject of our work is the neutral 
current QE neutrino scattering by nuclei, in what follows 
we would like to note various theoretical models that have 
been used also for the description of CC processes in 
connection to their applications to analyses of the NC 
processes. 

The analyses of the CCQE MiniBooNE data have 
raised questions on the capabilities of the various models 
to account for the different contributions to the neutrino- 
nucleus scattering cross sections. The RFG model, in 
which the shell structure and the nucleon correlations 
are neglected, gives results for the CCQE neutrino scat¬ 
tering that underestimate the data. The accordance with 
the data is achieved by increasing the world-average ax¬ 
ial mass M a {M a = 1.032 GeV) to M A = 1.35 GeV. 
An enhancement of the world-average axial mass is re¬ 
quired also in other models based on the impulse ap¬ 
proximation (IA) {e.g. P, [49M56I ] 4. This is an indica¬ 
tion that models based on the IA may lack important 
contributions to the processes of neutrino-nucleus scat¬ 
tering. In approaches beyond the I A, ingredients such as 
two-particle two-hole (2p — 2 h) contributions have been 
included. In the works [5711581 ] an approach based on the 


random-phase approximation, improved by considering 
relativistic corrections [HU led to a good agreement with 
the MiniBooNE data for both CCQE and neutral cur¬ 
rent quasielastic (NCQE) scattering, including the dou¬ 
ble differential CCQE cross section. It was pointed out 
in [53} that the multinucleon contribution may effectively 
be accounted for by increasing the value of the axial 
mass. As shown in [54 159l - l6l| . the RMF approach gives 
a good description of the shape of the double differen¬ 
tial cross section from the MiniBooNE experiment but 
fails to reproduce its normalization. It has been noted 
in [HU that meson exchange currents (MEC) could reduce 
the discrepancy. The calculations within the relativistic 
Green’s function (RGF) model [H, iH, [63j have provided 
a good description of the total CCQE and of the (dou¬ 
ble) differential (CCQE) NCQE MiniBooNE cross sec¬ 
tions [H, [Hi, EH, Ml, Ej]- The larger RGF cross sections 
can be attributed to the overall effect of inelastic chan¬ 
nels, such as, for instance, rescattering of the knocked- 
out nucleon, multinucleon processes, non-nucleonic A- 
excitation, that are recovered in the model by the use of 
a complex relativistic optical potential to describe FSI 
and that are not included in other models based on the 
IA. 

The SuperScaling Approach (SuSA) previoulsy dis¬ 
cussed has been used for analyses of neutrino-nucleus 
processes in [26|, H3, M, E3 Iii>- i66! ■ The results of SuSA 
underestimate by 20-30% the MiniBoonE data. An up¬ 
dated version of the model (SuSAv2), which incorporates 
different RMF effects in the longitudinal and transverse 
channels as well as isospin dependence, yields a milder 
disagreement (10-15%) (66|. The account for MEC in¬ 
creases significantly the results for the cross sections re¬ 
ducing the discrepancy with the data in the case of an¬ 
tineutrino but not so much in the neutrino case. 

Multinucleon effects on CC neutrino-nucleus scat¬ 
tering have also been investigated by the use of the 
Giessen Boltzmann-Uehling-Uhlenback event generator 
in Ref. [HU, where the energy spectra of the knockout 
nucleons are given in detail. 

In parallel, it is also interesting to compare with the 
older NCQE data from the Brookhaven BNL E734 exper¬ 
iment [44l ], corresponding to neutrino kinematics similar 
to MiniBooNE. We should emphasize that the NC cross 
sections depend also on the strangeness content of the nu¬ 
cleon, particularly, through the axial form factor. Hence 
a good control of nuclear effects is also required if one 
wants to extract information on the strange form factors 
of the nucleon from NC data. 

In [Hi} the differential and total cross sections for en¬ 
ergies ranging from a few hundreds of MeV to 100 GeV 
have been obtained and compared with the data from 
the BNL E734, MiniBooNKand NOMAD (see Ref. [H3). 
It has been concluded in [47| that the nuclear effects in 
NCQE and CCQE scattering seem to be very similar, 
though, according to the authors, the combined analyses 
of the CCQE and NCQE data does not seem to support 
the contribution of multinucleon final states being large 
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enough to explain the normalization of the MiniBooNE 
cross sections. It should be mentioned that in [53 an 
effective value of Ma = 1-23 GeV has been used. 

The sensitivity to FSI of NCQE (anti)neutrino scatter¬ 
ing cross sections has been investigated in @1i where the 
RGF cross sections calculated with different parametriza- 
tions for the phenomenological relativistic optical poten¬ 
tial are compared with the MiniBooNE data. The RGF 
results obtained with the standard value of the axial mass 
describe well the NCQE (anti)neutrino scattering data. 
It was pointed out, however, that the application of the 
RGF model to the semi-inclusive NCQE scattering can 
also include contributions of channels which are present 
in the inclusive but not in the semi-inclusive reaction. 

More theoretical approaches were presented in 
Ref. [5(|, where an axial mass Ma = 1-28 ± 0.05 GeV 
and a strangeness As = 0.11 ± 0.36 was extracted, and 
in Ref. |69fl within the transverse enhancement (TE) 
model |70j, 1711 1 , where a good agreement between the the¬ 
oretical results and the NCQE data was obtained, that 
the authors interpret as a very sizable contribution of 
2p — 2h and 3p — 3h processes. Finally, in Ref. [72j an ap¬ 
proach based on a realistic spectral function was applied 
to calculate the neutron-knockout (through neutral cur¬ 
rent scattering) cross sections in the kinematical regime 
of atmospheric-neutrino interactions in a broad energy 
range extending to 10 GeV. 

The main aim of this paper is to analyze neutral cur¬ 
rent QE (anti)neutrino scattering cross sections on a 12 C 
target using a realistic spectral function S(p,£) that gives 
a scaling function in accordance with the (e, e') scatter¬ 
ing data. In our previous work [T.ij this approach was 
applied to calculate CCQE (anti)neutrino scattering on 
12 C. The spectral function accounts for the NN corre¬ 
lations by using natural orbitals (NOs) from the Jas- 
trow correlation method and has a realistic energy de¬ 
pendence. In the calculations the standard value of the 
axial mass Ma = 1.032 GeV is used. As in [73], in the 
present work the role of FSI on the scaling functions and 
on the cross sections is taken into account. A compar¬ 
ison of the calculations with the empirical data of the 
MiniBooNE and BNL experiments is performed. 

The theoretical scheme of the work is given in Sec. El 
where the method to obtain a realistic spectral func¬ 
tion as well as the main relationships concerning the 
NCQE (anti)neutrino-nucleus reaction cross sections are 
presented. The results of the calculations and the discus¬ 
sion are given in Sec. lIIIl while the summary of the work 
and the conclusions are included in Sec. hyi 


II. THEORETICAL SCHEME 

The general formalism for NC (anti)neutrino scattering 
in the QE regime has been introduced in many previous 
works [3C)j . ]33t . IzUzB- Here we summarize briefly those 
aspects which are of more relevance for the later discus¬ 
sion of the results and of their comparison with Mini¬ 


BooNE and BNL data. We consider the semi-leptonic 
quasi-free scattering from nuclei in Born approximation, 
assuming that the inclusive cross sections are well repre¬ 
sented by the sum of the integrated semi-inclusive proton 
and neutron emission cross sections [30 ]. The kinematics 
for semi-leptonic nucleon knockout reactions in the one- 
boson-exchange approximation is presented in Fig. [I] 

A lepton with 4-momentum K^ = (e, k) scatters 
to another lepton with 4-momentum K ,fl = (e',k'), 
exchanging a vector boson with 4-momentum Q^ = 
K The lepton energies are e = V m 2 + k 2 and 
e' = y/m' 2 + k' 2 , where the masses of the initial and fi¬ 
nal lepton to and m' are assumed to be equal to zero 
for NC neutrino scattering. In the laboratory system 
the initial nucleus being in its ground state has a 4- 
momentum = ( M 0), while the final hadronic state 
corresponds to a proton or neutron with 4-momentum 
Pn=p or n = {En,Pn) and an unobserved residual nu¬ 
cleus with 4-momentum Pg = ( Eg,pg )■ Usually the 
missing momentum p = —p b and the excitation energy 

£ = Eb — Eg, with Eg = yj{Mg) 2 + p 2 are introduced, 
Mg being the ground-state mass of the daughter nucleus. 
Although in real situations, as is the case for MiniBooNE 
and BNL, there are usually no monochromatic beams and 
an integral over the allowed energies folded with the neu¬ 
trino flux must be performed, we assume the energy of 
the incident neutrino to be specified and also the outgo¬ 
ing nucleon energy En to be known. Finally, the angle 
dkp N between the incident neutrino and the ejected nu¬ 
cleon momentum is also given. 

Starting from the Feynman amplitude associated with 
the exclusive diagram of Fig. [TJ one can get inclusive 
cross sections by integrating over the undetected outgo¬ 
ing particles. 

In the case of QE electron- or CCQE neutrino¬ 
scattering, the outgoing lepton is detected and a sum 
over the outgoing nucleon variables is performed. Using 
the language introduced in this context in Ref. ;74|, we 
refer to these processes as “t-channel” reactions, since 




FIG. 1. Kinematics for semi-leptonic nucleon knockout reac¬ 
tions in the one-boson-exchange approximation. 
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the Mandelstam variable t = (K IJ — K'^) 2 is fixed. Then 
the (e, e') and CC ( vi,l ') reactions are governed by the 
same kinematics and the scaling formalism developed for 
the former can be trivially extended to the latter. 

In the case of NC neutrino scattering only the outgoing 
nucleon can be experimentally detected, while the unob¬ 
served outgoing neutrino is integrated over. This is re¬ 
ferred to as a “u-channel” process, where the Mandelstam 
variable u = (K^ — P At ) 2 is fixed. Then the kinematics is 
not the same as in the (e, e') case and, in particular, the 
two inclusive cross sections involve an integration over a 
slightly different region in the missing energy and mo¬ 
mentum plane. As a consequence it is not evident that 
the scaling arguments can be applied to NC scattering. 
However, in [30[ the influence of a non-constant in the 
derivation of scaling in the NC case was thoroughly inves¬ 
tigated within the general framework of the RFG model 
concluding that the scaling ideas still work properly for 
NC neutrino-nucleus processes. That study was extended 
in [zzllzi making use of the RMF approach. These re¬ 
sults showed that scaling of the second kind, i.e., inde¬ 
pendence of the nuclear target, works extremely well. On 
the contrary, scaling of first kind (independence on the 
transfer momentum) depends on the specific kincmatical 
situation considered. In general, first-kind scaling seems 
to be well respected when the angle of the ejected nucleon 
is larger than roughly 50°. This is the region where the 
cross section integrated over angles reaches larger values. 
Therefore, first-kind scaling is expected to work properly 
at MiniBooNE and BNL. Indeed, for the RMF and the 
particular kinematics involved in the experiments ana¬ 
lyzed in this work (MiniBooNE and BNL), we have veri¬ 
fied that the calculation of the NC cross sections based on 
u-scaling (as indicated in Eqs. m and (II 111 plus Ff) in 
the text below) gives rise to results very similar (within 
few percents) to those provided by the full calculation, 
i.e., without resorting to the scaling assumption. This is 
strictly true for Q 2 < 0.7 GeV 2 . For larger values of Q 2 , 
the u-scaling approximation begins to deviate from the 
full result, but by an amount significantly smaller than 
the uncertainty linked to the data error bands. 

The usual procedure for calculating the ( l,l'N ) cross 
section includes the Plane Wave Impulse Approximation 
(PWIA) and integrations over all unconstrained kine¬ 
matic variables. It is shown in [30j that the inclusive 
cross section in the u-channel can be written after some 
approximations in the following form: 


tion 


F («) 


1 1 
327T6 q' 



o 


x w^(p,p n)D V {Q 2 ) 2 


(3) 


is almost independent of (p, £) for constant (k,pN, @kp N )- 
In Eq. ([3]) and are the leptonic and s.n. hadronic 
tensor, respectively, and Dy(Q 2 ) is the vector boson 
propagator [13. In Eq. © — Pfc = (u/,q') 

is the four-momentum transferred from the initial lepton 
to the ejected nucleon and y' is the scaling variable nat¬ 
urally arising in the u-scattering kinematics, analogous 
to the usual ^/-scaling variable for f-scattering. The scal¬ 
ing function F(y') obtained within a given approach can 
be used to predict realistic NC cross sections. Assuming 
that the domains of integration T> u (in the u-channel) 
and T>t (in the f-channel) are the same or very similar, 
the results for the scaling function obtained in the case 
of inclusive electron scattering (where T>t works) can be 
used in the case of NC neutrino reactions. It is pointed 
out in [13 that V t and V u differ significantly only at large 
£ (also at large p , but there the semi-inclusive cross sec¬ 
tions are expected to be negligible). So, given that the 
semi-inclusive cross sections are dominated by their be¬ 
havior at low £ and low p, one expects the results of the 
integrations in the t- and u-channel to be very similar, 
and thus the scaling functions will be essentially the same 
in both cases. 

As noted in [30], if the s.n. cross section is smoothly 
varying within the (p, £) integration region, the differen¬ 
tial cross section in the RFG can be factorized as shown 
in Eq. m with the RFG scaling function: 


Frfg(V’RFg) = 2 kp l 1 - ^RFg) 0 (! ~ ^rfg) > ( 4 ) 

where the RFG it-channel '0-variable is defined: 


4>rfg = s 



1 + 


/ Urfg 2 

V mjv J 


1/2 



(5) 


and 


tun 

11 rfg = s —— 


X' \/ T' 2 p ' 2 + t' — k't' p' 


( 6 ) 


da 

dfl^dpN 




where 


F{y',q') = J pdp J ~ F{y'), 

v u 


( 1 ) 


( 2 ) 


provided the effective NC single nucleon (s.n.) cross sec- 


is the RFG y-scaling variable for the u-channel and cor¬ 
responds to the minimum momentum required for a nu¬ 
cleon to participate in the NC neutrino-nucleus scatter¬ 
ing. The dimensionless kinematic quantities in Eq. © 
are given by: n' = q' /2m,N, X' = u // 2 toat , t ' = k' 2 — A' 2 , 

p' = 1 ——(1 — m' 2 /m 2 N ). The sign s is 


s = sgn 


1 

t' 


X'\J t' 2 p ' 2 + t' 


k't' p 



(7) 















5 


The basic relationships used to calculate the s.n. cross 
sections are given in [30(. This concerns the leptonic and 
hadronic tensors and the response and structure func¬ 
tions. Here we summarize the basic expressions for the 
neutral weak nucleon form factors and the specific set of 
parameters considered in the calculations: 

• Weak mixing angle: sin 2 9w = 0.23122(15), [79]; 

• Isovector axial form factor of the nucleon at zero 
momentum transfer: g A = 1.2695, [7§|; 

• The weak form factors are given in the form 
[Eqs. (55) and (56), Ref. My 


S(p,S) [10’ 9 MeV 3 ] 



0 50 100 150 200 250 300 

p [MeV] 


G a {Q 2 ) = ^ =1 Gl =1 r 3 + V3#=°G£= 0 + , 

a = {E,M}, 

Ga(Q 2 ) = & =1 Gfn + e A =°G ( 2 ] + d 0) Gi s) , 

with (no radiative corrections included) 

£ y =1 = 2(1 - 2 sin 2 9 W ); = —4 sin 2 9 W ; 

{f o) = -i ; d- I = -2;S-“ = 0;{r = i. 

• The isoscalar and isovector form factors are (for 
more details, see ®): 

/^iT—1 s^ip /-in /-iT —0 (~iV s~in 

^E,M ~ ^E,M ^E,Mi ^E,M ~ ^E,M ^E,Mi 

Gf = g A /2 GUQ 2 ), G^ = «?i s) Gp(Q 2 ), 
with G^(Q 2 ) = (1 + \Q 2 \/M\)~ 2 . 

• Electric and magnetic strange form factors (dipole 
Q 2 -dependence has been used): 

G ( e\Q 2 ) = PstG^Q 2 ), G$(Q 2 ) = p s G v D {Q 2 ). 

• The electromagnetic (EM) form factors correspond 
to the so-called GKex (Gari-Kriimpelmann ex¬ 
tended, developed by E. L. Lomon) model Bp IsTif . 

• Unless otherwise specified, the following set of 
strange parameters has been considered (see 
Ref. [8(| for more details): 

H. = -0.020; p s = 0.59; g = 0. 

• Finally, the axial mass is taken to be the world- 
average value Ma = 1.032 GeV. 

In this work we use in the calculations the RFG, 
SuSA, harmonic oscillator (HO+FSI) and natural or¬ 
bitals (NO+FSI) scaling functions including final state 
interactions (see [73} and the text below). We present 
also results obtained in the RMF and RGF models, where 
the cross sections are calculated in a fully unfactorized 
approach which does not make use of the approxima¬ 
tions leading to Eq. ©■ All our results are used to ana¬ 
lyze NCQE (anti)neutrino cross sections on a CH 2 target 


FIG. 2. (Color online) The 12 C realistic spectral function 
S(p,£), which is constructed using natural orbits single¬ 
particle momentum distributions from the Jastrow correlation 
method and Lorentzian function for the energy dependence 
(see text). 


measured by the MiniBooNE collaboration [43|, |46| . We 
shall also compare with the BNL E734 experiment H |. 
studying vp and Vp NCQE interactions, where the target 
was composed in 79% of protons bound in carbon and 
aluminum and in 21% of free protons. 

In details, here we present briefly how the HO+FSI 
and NO+FSI scaling functions are obtained (see also 
Ref. [Zl): 

(i) The spectral function S(p,£) is constructed in the 
form 

S(p,£) = '^2‘2(2ji + l)ni(p)L Ti (£ ~ ( 8 ) 


(ii) The single-particle momentum distributions rq(p) 
are taken to be either corresponding to the HO 
single-particle wave functions or to the NO’s from 
the Jastrow correlation method, where central 
short-range NN correlations are included; 


(iii) The Lorentzian function 


Lr<{£-£i) 


i r,/2 
*{£- £i? + (iy2) 2 ’ 


(9) 


is used for the energy dependence of the spectral 
function with parameters Ti p = 6 MeV, Tis = 
20 MeV, which are fixed to the experimental widths 
of the Ip and Is states in 12 C. 

The realistic spectral function S(p,£) is presented 
in Fig. [2j where the two shells lp and Is are clearly 
visible. 


(iv) For a given momentum transfer q and energy of the 
initial electron e we calculate the electron-nucleus 
( 12 C) cross section using 


da t 

dcod\q\ 



dE d 3 p 


St(p,E) 

EpE p i 


x6(w + M-E-Ep.)L%HZ it 


( 10 ) 








6 


in which the spectral function S(p,£) [Eq. (0] is 
used, the index t denotes the nucleon isospin, L®“ 
and t are the leptonic and hadronic tensor, 
respectively (for details, see Ref. 0); 


(v) The corresponding scaling function F(q,uj) is cal¬ 
culated within the PWIA by 


F{q,u) 


[da/de'dd '] (ee ,) 
a eN (q,u-,p= \y\,£ = 0) 


( 11 ) 


and by multiplying it by kp the scaling function 
f(ip) is obtained, where kp is the Fermi momen¬ 
tum for a specific nucleus and lf eN is the azimuthal 
angle-averaged single-nucleon cross section; 


(vi) To account for FSI, the (5-function in Eq. m is 
replaced by 


S(u + M-E-Ep>) 


W/ 7 


W 2 + [u> + M — E — E p > — V] 2 


( 12 ) 


where the real ( V ) and i mag inary (W) parts of the 
OP are obtained in Ref. [84[ from the Dirac OP. 


III. RESULTS AND DISCUSSION 


In this section the theoretical predictions of the RFG, 
HO+FSI, NO+FSI, SUSA scaling functions are com¬ 
pared with the data measured by the MiniBooNE and 
BNL Collaborations. The comparison is performed also 
with the results of the RMF and RGF models, which 
are based on the same relativistic mean-field model for 
nuclear structure but on a different treatment of FSI. In 
the RMF model FSI are described by the same relativistic 
mean field potential describing the initial nucleon state; 
the description of FSI in the RGF is based on the use of 
a complex optical potential. Details of the RGF model 
can be found, for instance, in (gH, [86]. The results of 
the RMF and RGF models have been already and widely 
compared in [63| for the inclusive QE electron scatter¬ 
ing, in [55], [III for CCQE and in [61| for NCQE neutrino 
scattering. The RGF calculations presented in this work 
have been carried out with the recent democratic optical 
potential of [H| (RGF-DEM). 

The comparison between theory and experiment for 
the NCQE flux-averaged MiniBooNE (anti)neutrino 
cross section is presented in Fig. [3] Here we compare the 
predictions using the RFG, HO+FSI, NO+FSI, SUSA 
scaling functions, and RMF model with the data. As 
usual in NC reactions, in this work, the variable Q 2 is 
defined as Q 2 = 2 M^Tn, where Mn and Tn are the 
mass and kinetic energy of the outgoing nucleon, respec¬ 
tively. In order to compare with MiniBooNE we evaluate 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

Q 2 [GeV 2 ] 


FIG. 3. (Color online) NCQE neutrino [panel (a), vN —» 
uN] and antineutrino [panel (b), VN —► VN] flux-averaged 
differential cross section computed using the RFG, HO+FSI, 
NO+FSI, SUSA scaling functions, RGF and RMF models and 
compared with the MiniBooNE data [43 . [46]. The results 
correspond to the world-average axial mass Ma = 1.032 GeV 
and strangeness As = 0. The error bars do not account for 
the normalization uncertainty of 18.1% (19.5%) in the v(y) 
case. 


the following differential cross section per nucleon: 


d<J v (p)N^v{p)N 1 


dQ 


r< d<Jy(v)p^v(v ) Pl H 

2 n ) 


dQ 2 
dcr 1, (j7) p—t 1/(77} p^c 


c(<3 * dQ 2 

■l0^ n , ci Q 2 ) ^^S^ , (13) 


where C' y(77)PiH (Q 2 ), C^p.cCQ 2 ), and C^w.cf Q 2 ) are 
the efficiency correction functions, given in Refs. [43l. |46I| . 
of three different processes: the (anti)neutrino scattering 
off free protons in the hydrogen atom, the bound protons 
in the carbon atom and the bound neutrons in the carbon 
atom, dcv^p^p^u^p^pii/dQ , dcy l/ (jd)p—ru(j7)p.c./d(J , and 
da v {T 7 ) n ^nd(v)n,c/dQ 2 are the theoretical i/(F)NCQE cross 
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FIG. 4. (Color online) Differential cross section da/dQ 2 for 
I'NCQE scattering per bound nucleon in the 12 C nucleus cal¬ 
culated within RMF model and multiplied by the MiniBooNE 
flux. The results are given at several fixed values of four- 
momentum transfer squared (Q 2 ). 



FIG. 5. (Color online) Ratio of the antineutrino to neutrino 
NCQE scattering cross section computed using the RFG, 
HO+FSI, NO+FSI, SUSA, RGF, and RMF approaches and 
compared with the MiniBooNE data [46). 


sections on free protons (per free proton), on bound pro¬ 
tons (per bound proton), and on bound neutrons (per 
bound neutron). As can be seen from Fig. [3] the the¬ 
oretical results corresponding to all models except the 
RGF-DEM underestimate the neutrino data in the re¬ 
gion between 0.1 < Q 2 < 0.7 GeV 2 , while all theories are 
within the error bars for higher Q 2 . On the other hand 
the same models underestimate the antineutrino data at 
high Q 2 . This is clearly seen in the insets of Fig. [3l where 
the cross sections are represented in logarithmic scale. 
The RGF-DEM results are larger than the results of the 
other models and in generally good agreement with the 
data over the entire Q 2 region considered in the figure. 
The enhancement of the RGF cross sections is due to the 
contribution of final-state channels that are recovered by 
the imaginary part of the optical potential and that are 


not included in the other models. 

In order to better understand the behaviour of our re¬ 
sults with Q 2 , we present in Fig. 0] the differential cross 
sections da/dQ 2 for i^NCQE scattering (per bound nu¬ 
cleon in the 12 C nucleus) calculated within the RMF 
model and multiplied by the MiniBooNE flux, at sev¬ 
eral fixed values of Q 2 . <h t ot is the total integrated 
flux factor for the MiniBooNE experiment: 

$tot = J §(£v)d£ v . (14) 

As can be seen from Fig.[I]for small values of Q 2 the con¬ 
tribution to the cross sections is accumulated from a large 
neutrino energy range (for example, at Q 2 = 0.19 GeV 2 
the contributions to the cross section come from neu¬ 
trino energies from 0.2 to 2.5 GeV), whereas for higher 
Q 2 values, the energy range which gives contribution to 
the cross section becomes smaller (for example, at Q 2 = 
1.69 GeV 2 contributions to the cross section come from 
neutrino energies from 1.1 to 2.5 GeV). This suggests that 
the discrepancy between the RMF result and the exper¬ 
imental data at low Q 2 can be ascribed to the failure of 
the impulse approximation for low (anti)neutrino ener¬ 
gies. For the other models (leaving apart the RGF) the 
situation is different for v and v scattering. To better 
illustrate this difference we show in Fig. [5] our predic¬ 
tions for the ratio of the antineutrino to neutrino NCQE 
scattering cross sections. The results are compared with 
the recently published MiniBooNE data 0] . Since both 
measurements were made in the same beamline and with 
the same detector, it is expected that the bin-to-bin ratio 
of the two cross section measurements would cancel the 
common systematic errors. As can be seen from Fig. [5] 
the results of all models are in agreement with the data 
up to Q 2 < 0.8 GeV 2 (actually the RMF and RGF results 
are in good agreement up to Q 2 ~ 1 GeV 2 ), whereas for 
bigger Q 2 all models underpredict the data. This is so 
because the theoretical antineutrino NCQE cross section 
underpredicts the data at high Q 2 (see Fig. n. It has to 
be noted that the results using different scaling functions 
almost coincide, since the cross sections are proportional 
to s.n. cross sections multiplied by the scaling function 
[see Eq. ©] and the ratio is proportional to the ratio 
of the s.n. antineutrino to s.n. neutrino cross sections. 
The updated version of the SuSA model, “SuSAv2”, gives 
an interesting possibility in the studies of the NC cross 
sections and particularly of the mentioned ratio. This 
is related to the use of different parametrizations of the 
transverse and longitudinal scaling functions in SuSAv2. 
Work along this line is in progress. 

For completeness we also present in Fig. [6] the spectra 
corresponding to the numerator and denominator enter¬ 
ing the ratio between ^-scattering from proton and nu¬ 
cleon (proton plus neutron) (left panels) and the ratio 
computed by dividing of the two samples (right panels) 
within the various models and compared with the Mini¬ 
BooNE data . The numerator and denominator data 
(left panels of Fig. [6]) are taken from 0], where the data 
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FIG. 6. (Color online) The RFG, HO+FSI, NO+FSI, SUSA, RGF, and RMF predictions, after the folding procedure, compared 
with the histograms of the numerator (top-left panel) and denominator (bottom-left panel) of the ratio. The error bars in the 
left panel represent only the statistical uncertainty, computed as the square root of the event number. The corresponding ratio 
is shown in the right panel of the figure. The axial mass and strangeness have been assumed to be the standard axial mass 
value and zero strangeness. Data are taken from f43j, [&!, [90|. 




FIG. 7. (Color online) NCQE flux-averaged cross section: up —> up [panel (a)] and up —y up [panel (b)[ compared with the 
BNL E734 experimental data [44j]. Our results are evaluated using the RFG, HO+FSI, NO+FSI, SUSA scaling functions, and 
RMF model with the standard value of the axial-vector dipole mass Ma = 1.032 GeV 2 and strangeness As = 0. The error 
bars do not include the normalization uncertainty of 11.2% (10.4%) in the u(u) case. 


are reported without the corresponding errors (so, in the 
figure only statistical errors are included). In the calcula¬ 
tions the axial mass and strangeness have been assumed 
to be the standard axial mass and zero strangeness. We 
note that the dispersion between the models tends to 
cancel when this ratio is considered. This result clearly 


shows that the proton/nucleon ratio is very insensitive 
to nuclear model effects and to final state interactions, 
and hence, it may provide information that improves our 
present knowledge on the electroweak nucleon structure, 
in particular, the nucleon’s strangeness. In particular, for 
kinetic energies of the outgoing nucleon X)v > 350 MeV 
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our models give results which are in good agreement with 
the experimental data (left panels of Fig. [6]), while for the 
ratio our theoretical results are within the error bars at 
all kinematics (right panel of Fig. IS). 

We now compare the results obtained with our models 
with the BNL E734 experimental data. The mean value 
of neutrino (antineutrino) energy is 1.3 GeV (1.2 GeV) 
for BNL experiment, while for MiniBooNE experiment 
it is 788 MeV (665 MeV). In Fig. [7] the differential cross 
sections evaluated using the RFG, HO+FSI, NO+FSI, 
SUSA scaling functions, and RMF model are com¬ 
pared with NCQE vp —> vp [panel (a)] and vp —> vp 
[panel (b)] BNL E734 experimental data. The BNL 
E734 experiment was performed using a 170-metric- 
ton high-resolution target-detector in a horn-focused 
(anti)neutrino at the Brookhaven National Laboratory. 
The cross sections results show a behaviour similar to 
those of the MiniBooNE experiment. The latter (using 
the Cherenkov detector filled with mineral oil (CH 2 )) 
is sensitive to both v(y)p and v(y)n NCQE scatter¬ 
ing SLIM- It has been known for some time (see, 


e.g., ZS~Zfil) that the As-dependence of the NCQE 
neutrino-nucleon cross section is very mild. This results 
from a cancellation between the effect of As on the proton 
and neutron contributions, which are affected differently 
by the axial strangeness: by changing As from zero to 
a negative value the proton cross section gets enhanced 
while the neutron one is reduced, so that the net effect 
on the total cross section is very small. NCQE v(v)p dif¬ 
ferential cross sections were measured in the BNL E734 
experiment, which are sensitive to the values of As (there 
is not a cancellation effect). The BNL E734 experimen¬ 
tal data can be reproduced within our models in principle 
by the fit of the axial strangeness without change of the 
axial mass value. 

Here we would like to mention that, first, our calcula¬ 
tions using NO and HO single-particle wave functions in 
rii(p) in the spectral function Eq. (|5|) with FSI and with¬ 
out FSI show that the inclusion of FSI effects leads to a 
small change (a depletion) of the cross sections. Second, 
the results for the cross sections obtained using realis¬ 
tic spectral function S(p,E) with single-particle momen¬ 
tum distributions rii(p) (see Eq. ([8])) that include Jastrow 
short-range NN correlations (accounted for in the NO’s) 
can be compared in Figs. [3] [5j [6] and [7] with those when 
NN correlation are not included (RFG and HO). It can 
be seen that, similarly to the case of CCQE neutrino 
scattering (see Ref. [7J]), the differences between results 
in correlated and non-correlated approaches are small, 
thus showing that the process is not too sensitive to the 
specific treatment of the bound state. 


IV. CONCLUSIONS 

This work complements previous studies focused on 
charged-current quasielastic (CCQE) neutrino-nucleus 
scattering processes that were performed making use of a 


realistic spectral function. Here we extend this analysis 
to the case of neutral-current (NC) neutrino processes, 
and compare our theoretical predictions with data mea¬ 
sured by MiniBooNe and BNL Collaborations. Contrary 
to CCQE reactions, where the final lepton is detected, in 
the NC case one has no information on the energy and 
momentum of the ejected neutrino. Hence the trans¬ 
ferred four-momentum cannot be determined. This, as 
already discussed in some previous works, makes the de¬ 
scription of the reaction mechanism not so clear and some 
caution should be also drawn on the “validity” of scaling 
arguments when applied to NC. However, our previous 
studies give us confidence in the reliability of our cal¬ 
culations and their application to describe the present 
experimental data measured at different facilities. 

The main objective of this work is centered in the use 
of a realistic spectral function, that accounts for short- 
range NN correlations, and has also a realistic energy 
dependence. This function gives a scaling function in ac¬ 
cordance with electron scattering data and it can be used 
for a wide range of neutrino energies. Therefore, the use 
of this spectral function to describe the general reaction 
mechanism involved in NC neutrino-nucleus scattering 
processes can provide very valuable information that can 
be confronted with results obtained with other theoret¬ 
ical approaches. In this sense, we compare our spectral 
function-based predictions with the results provided by 
the SuSA, RMF and RGF models largely used by us in 
the past. The discrepancies found can help disentangling 
effects directly linked to particular ingredients in the pro¬ 
cess: final state interactions, nucleon correlations, effects 
beyond the impulse approximation, etc. 

The predictions of our model agree in general with pre¬ 
vious results although with some peculiarities that should 
be emphasized. Our calculations showed that the inclu¬ 
sion of FSI effects in the spectral function-based calcula¬ 
tions leads to a slight depletion of the cross section being 
in close agreement with the RFG prediction. The in¬ 
clusion of FSI effects in the RGF model leads to larger 
cross sections, in good agreement with the data. On the 
contrary, SuSA and, in particular, RMF approaches lead 
to significantly smaller differential cross sections at low 
values of |Q 2 | (< 0.6-0.8 GeV 2 ) departing also from the 
data. This behavior can be seen for the two experiments 
considered in the work: MiniBooNe and BNL. Another 
point of relevance when comparing the different models 
is the softer Q 2 dependence (with a smaller slope) shown 
by the RMF cross section. Whereas it is clearly below 
the other curves at low Q 2 (up to ~ 0.5-0.6 GeV 2 ), it 
crosses them providing the largest contribution at higher 
Q 2 . This result can be taken as an indication of the 
particular sensitivity of NC processes to the specific de¬ 
scription of FSI effects. It may be also connected with 
the increasing tail in the scaling function provided by the 
RMF model at larger Q-values. This is a consequence 
of the enhancement of the lower components in the rela¬ 
tivistic nucleon wave functions, particularly, for the final 
state. 
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All our calculations are based on the impulse approxi¬ 
mation, i.e., they do not include effects beyond the one- 
body approach, for example 2 p — 2h contributions in¬ 
duced by meson exchange currents. These ingredients 
have been shown to be very important in the analysis of 
neutrino-nucleus scattering processes. In particular, they 
produce a significant enhancement in the cross section at 
low-to-moderate values of the transferred four momen¬ 
tum. This is consistent with our predictions that clearly 
underestimate data for such kinematical regions. On the 
contrary, the accordance improves at higher Q 2 . This 
is strictly true for neutrinos at MiniBooNE. In the case 
of antineutrinos, MiniBooNE data at Q 2 > 1—1.2 GeV 2 
are higher than theoretical predictions, the RMF results 
being closer to the experiment. This behavior leads also 
to the significant discrepancy observed for the antineu¬ 
trino/neutrino ratio (Fig. [5J. More studies are needed 
in order to understand these differences at medium-large 
<5 2 -values. This could be related to a different role played 
by 2 p— 2 h contributions and MEC for neutrinos and an- 
tineutrinos. Work is in progress to evaluate the impact 
of 2 p — 2 h excitations on these results. 
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